Informative priors and the analogy between quantum and classical heat engines 
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When incomplete information about the control parameters is quantified as a prior distribution, 
a subtle connection emerges between quantum heat engines and their classical analogs. We study 
the quantum model where the uncertain parameters are the intrinsic energy scales and compare 
with the classical models where the intermediate temperature is the uncertain parameter. The prior 
distribution quantifying the incomplete information has the form n(x) oc 1/x in both the quantum 
and the classical models. The expected efficiency calculated in near-equilibrium limit approaches 
the value of one third of Carnot efficiency. 
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I. INTRODUCTION 

Quantum heat engines (QHEs) are novel tools to study 
the underlying thermodynamic properties of quantum 
systems The working substance in a QHE is a 

few-level quantum system and can show interesting fea- 
tures like quantum correlation, coherence and so on. So 
QHEs may show unexpected behavior [l|, 0, H[ which is 
not possible in the classical models of heat engines. But 
all these models are consistent with the second law of 
thermodynamics. Recent studies @, [|[ showed that the 
expected behavior of certain models of QHE exhibit clas- 
sical thermodynamic features which points out interest- 
ing and novel connection between information and ther- 
modynamics. In these models, the uncertain parameters 
are treated as in Bayesian approach. 

In Bayesian approach to probability theory, prior dis- 
tribution [To| or known simply as a prior, quantifies 
the prior knowledge about the uncertain parameter(s). 
Usually, there is some information available even about 
the uncertain parameters, e.g. from the nature of the pa- 
rameters or from the physics of the problem. The prior 
which makes use of this knowledge is also addessed as an 
informative prior. So the right choice of the prior plays 
an important role in this approach. 

In this paper, we discuss a quantum and a classical 
model of heat engine and estimate their performance. In 
the quantum model we consider a pair of two level sys- 
tems with energy level spacings a\ and et2. Reservoirs 
associated with the respective systems are at tempera- 
tures T\ and T2. These intrinsic energy level spacings can 
be controlled externally e.g. through external magnetic 
field. In the classical model, the pair of two level systems 
is replaced by a pair of classical ideal gas systems. 

In the quantum model, the unknown parameters are 
the energy level spacings of the two level systems. But in 
the classical case, the uncertain parameter is the interme- 
diate temperature. To assign the prior, we invoke differ- 



ent observers who satisfy a consistency criterion and thus 
arrive at the prior for the unknown parameter. These 
prior distributions are used to estimate the expected be- 
havior of thermodynamic quantities. Finally we compare 
the estimated values of the physical quantities obtained 
from the quantum and the classical models. The main 
objective of this paper is to show the equivalence of the 
expected behavior of quantum and classical models under 
certain conditions. Interestingly the expected efficiencies 
are also related to the efficiencies at optimal performance 
for certain finite time models of Brownian heat engines. 

The paper is organised as follows. In Section II, we 
present the quantum model for heat engine and sum- 
marise its main features. In Section III, we discuss the 
prior chosen in the quantum model based on the initial 
incomplete information about the internal energy scales 
of the working medium. In Section IV, we apply the 
prior so derived to estimate internal energies of the sys- 
tem. Further, we highlight a specific asymptotic limit in 
which the expected behavior becomes especially simple. 
Section V, we introduce the classical model with interme- 
diate temperature as the uncertain parameter. The final 
Section VI is devoted to conclusions and future outlook. 



II. THE QUANTUM MODEL 

Consider a pair of two-level systems labeled R and S, 
with hamiltonians Hr and H$ having energy eigenvalues 
(0,ai) and (0,et 2 ), respectively. The hamiltonian of the 
composite system is given by H = Hr ® I + i" ® Hg. 
The initial state is p; n i = pr ® ps, where pr and ps 
are thermal states corresponding to temperatures T\ and 
Ti (< Xi), respectively. Let (ri,ra) and (si,S2) be the 
occupation probabilities of each system, where 

1 1 m 

1 (l + e -ai/Tx)' 51 (l + e-to/Ti)' W 

with f2 = (1 — ri) and S2 = (1 — Si). We have set 
Boltzmann's constant fee = 1. The initial mean energy 
of each system is 
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where i — 1,2 denote system R and S respectively. 
Based on quantum thermodynamics [TTl - fl3j , the process 
of maximum work extraction is identified as a quantum 
unitary process on the thermally isolated composite sys- 
tem. It preserves not just the magnitude of Von Neu- 
mann entropy of the composite system, but also all eigen- 
values of its density matrix. It has been shown in these 
works that for a\ > a 2l such a process minimises the final 
energy if the final state is given by p& n = ps ^ PR- Ef- 
fectively, it means that in the final state the two systems 
swap between themselves their initial probability distri- 
butions. The final energy of each system at the end of 
work extracting transformation is 
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where i =/= j. The average heat absorbed from system R 
is defined as Q\ = Tr[(pfl — ps)Hr], is given by 



Qi 
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(4) 



Similarly, the average heat released to system S is defined 
as Q 2 = Tr[(p s - p R )H s ], is given by 



Q2 



-0,2 
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(l + e «i/ri) (i + e a 2 /T 2 ) 



(5) 



The average work done in one cycle is W = Q\ + Q 2 . To 
complete the cycle, the two systems are brought again 
in thermal contact with their respective reservoirs. The 
operation of the machine as a heat engine implies W > 
and Q\ > 0, which is satisfied if 



ai(r 2 /Ti) < a 2 < 01. 



III. PRIOR DISTRIBUTION 



(6) 



engine ? We have suggested a subjective or Bayesian 
approach to address this question @,@|. This implies that 
an uncertain parameter is assigned a prior distribution, 
which quantifies our preliminary expectation about the 
parameter to take a certain value. Thus the prior should 
be assigned by taking into account any prior information 
we possess about the parameters. For example, if at is 
specified, then the prior distribution for a 2l 7r(a 2 |ai) is 
conditioned on the specified value of a\, and is defined in 
the range [ai6, a{\, where 9 = T 2 /Ti, because we know 
the set-up works like an engine if we implement Eq. ([5]) . 
We denote the prior distribution function for our problem 
by II(ai, a 2 ). To assign the prior, it seems convenient to 
involve two observers A and B, who wish to assign priors 
for ai and a 2 . Based on the derivation given in [8j, we 
find that the prior for each parameter is 



n(oi) = 



In 



(7) 



and the joint prior for the system acting as an engine, is 
given by 



n(ai,a 2 ) 



1 



1 



ln (i) ln fc) aia2 ' 
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IV. EXPECTED VALUES OF QUANTITIES 

In this section, we use the priors assigned above, to find 
expected values for various physical quantities related to 
the engine. The expected value of any physical quantity 
X which may be function of a\ and a 2 , is defined as 
follows: 



X = 



X II(oi, a 2 ) daida 2 . 



(9) 



Now consider a situation in which the temperatures of 
the reservoirs are given a priori such that T\ > T 2 , but 
the exact values of parameters a\ and a 2 are uncertain. 
The prior information about these parameters may be 
summarised as follows: 

• a\ and a 2 represent the same physical quantity, i.e. 
the level spacing for system R and S respectively, 
and so they can only take positive real values. 

• If the set-up of R + S has to work as an engine, 
then criterion in Eq. (J6j) must hold, whereby if one 
parameter is specified, then the range of the other 
parameter is constrained. 

Apart from the above conditions, we assume to have no 
information about a± and a 2 . The question we address 
in the following is: What can we then infer about the 
expected behaviour of physical quantities for this heat 



These expected values reflect the estimates by an ob- 
server who assigns the priors. In principle, there are two 
ways to calculate the expected value of some quantity 
which depends, in general, on the method used. The 
method by which the observer A applies the joint prior 
is based on the definition 

n{a 1 ,a 2 )=n{a 2 \a 1 )n{a 1 ), (10) 

i.e. the prior for a\ is assigned first, followed by II(a2|ai) 
which is the conditional distribution of a 2 for a given 
value of a\. On the other hand, observer B applies 

n(a x ,a 2 ) =n(o 1 |o2)n(o 2 ) J (11) 

whereby the prior for a 2 is assigned first and n(ai|a 2 ) 
represents the conditional distribution of a\ for a given 
value of a 2 . 
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A. Internal energy 



B. Asymptotic Limit 



We calculate the expected values of internal energies 
for systems R and S. These values can then be used to 
find the expected work per cycle, heat exchanged and so 
on. 

(i) Initial state: For a given ai, the internal energy 
E^l is given by Eq. J2J). The expected initial energy is 
denned as 



^ini - 



(i) 

where i = 1,2. Note that E- ln [ depends only on a^, so we 
need to average over the prior for m only. Using Eqs. ([2]) 
and ([7]), we obtain 



(12) 
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1 + e a -"/ T » 
1 + e «w/Ti 
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(ii) Final state : In this case, the internal energy of 
R as well as S, is function of both a\ and a 2 (see ©) 
and so the expected values are obtained by averaging 
over the joint prior, II(ai, a 2 ). For instance, the expected 
final energy of system S (denoted by superscript (2)) as 
calculated by A is 

E^{A) = J J ^4n ) n(a 2 |ai)n(ai)da 1 da 2 



K 



ai 

dai I da 2 , 



K(l-d) 



Tiln 



(14) 



where K = [In (1/9) In (a max /a m ; n )] . Similarly as cal- 
culated by B, we have 



4? (B) 



J J E^U(a 1 \a 2 )U(a 2 )da 1 da 2 



K 



da-i 



a.2/8 



dai 



(1 + e Q i/ T i)ai 



(15) 



which cannot be solved analytically. 

Now in general, the expected final energies of S, as 
given by Eqs. (fT4|) and (fTo*]) according to A and B, re- 
spectively, are not equal. One would expect that if the 
state of knowledge of A and B is similar, then they should 
expect the same value for a given quantity. (Similar fea- 
ture is also observed in the expressions for expected final 
energy of system R.) 



As remarked above, observers A and B should arrive 
at similar estimates for physical quantities using their 



respective priors, 
and a max » Ti. 



This happens in the limit, 



<< To 



Then, Eq. (|13[) is approximated as 
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(16) 



The ratio (a ma x/amin) in the above may be large in mag- 
nitude, but is assumed to be finite. 

Similarly, considering the final energies, it is remark- 
able to note that in this limit, not only the expected 
energy of a system (R or S) calculated by either of the 
methods (A or B), is the same but also its value for sys- 
tem R or S is also equal. Thus from Eqs. (fT4"|) and (fl5|) . 
we have (omitting the observer index) 



E 



In 2 (l-6»)Ti 
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Further insight may be obtained if we estimate the final 
temperatures of systems R and S, after the work extrac- 
tion process. Note that if values of both a\ and a 2 are 
specified, the temperatures (T') of the two systems after 
work extraction, are given by [13j 



T;=T 2 ^, and T'=TA 



a 2 



«i 



(18) 



In general, the two final temperatures are different from 
each other. However following the subjective approach, 
when we look at the expected values of the final temper- 
atures as calculated by A or B, we find 



T 1= T 2 = Ti 



(1-0) 
m(l/6>)' 



(19) 



It is interesting to find that the assignment of the prior is 
such that the two systems are expected to finally arrive 
at a common temperature. Going back to Eqs. (TT6l) and 
(IT71) for the energies, we see that they satisfy a simple 

relation E^ cx Ti and E^ n oc t\. This is analogous to 
the thermodynamic behavior of a classical ideal gas. 

Next, the expected values of the heat exchanged be- 
tween system i and the corresponding reservoir is given 

by — E^ — Eft n . Q { > (Q { < 0) represents heat 
absorbed (released) by the system. Then the expressions 
for the heat exchanged with the reservoirs in the said 
limit are as follows: 

'" 2 ' 1+ (Izi>K (20) 
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ln2 ( (l-0)\ 

1 + K^4- I T 2 . (21) 
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Now the expected work per cycle is defined as: W — Q 1 + 
Q 2 - Thus the efficiency may be defined as r\ = 1 + Q 2 /Qi- 
Explicitly, using Eqs. (|20| and (f2Tj) we get 



»7 



1 



» In 6>+ (1 



In 6>+ (1 



(22) 



This is the efficiency at which the engine is expected to 
operate for a given 9. The above value is function only 
of the ratio of the reservoir temperatures. 

We note that the constant of proportionality in Eqs. 
(fl6|) and (JTTJ) , which is ln2 • (ln(a m ax/a m in))~ , can be 
related with heat capacity. The expected value of initial 
heat capacity of system i, defined as 



Ci = 



(23) 



where we know for a two level system, the canon- 
ical heat capacity at constant volume is Ci = 
(Oi/iy^exp (a i /T i ))/(l + exp {ajT,,)) 2 . 

In particular, for the asymptotic limit, the leading term 
yields 



Ci = C ! 



In 2 



In 



(ST 



(24) 



This value is independent of temperature of the system 
and thus indicates an analogy with a constant heat ca- 
pacity thermodynamic system. 

Thus the requirement of consistency between the re- 
sults of A and B implies, in an asymptotic limit, that 
the behavior expected from minimal prior information is 
the one which shows simple thermodynamic features such 
as constant heat capacity and equality of subsystem tem- 
peratures upon maximum work extraction. 



V. CLASSICAL MODEL 

In this section, we discuss a classical model of the heat 
engine, within the subjective approach. Consider two 
thermodynamic systems at initial temperatures 7\ and 
T% (< 2i), such that the heat capacity C can be assumed 
to be constant i.e. the systems behave like classical ideal 
gases. Further consider the maximum work extraction 
by coupling the two systems to a reversible work source. 
This process preserves the total entropy, i.e. AS = 0. 
Let us assume that at some stage, the temperatures of 
the systems are T a and T b , respectively. Entropy conser- 
vation leads to 



n = 



TiT 2 



(25) 



The above equation relates Tb and T a , implying that 
given a value of one of them, the value of the other is 
fixed. The work W extracted from the system is given 
by the decrease in internal energy, 



W = C(T 1 +T 2 -T a -T b ). 



(26) 



Similarly, the heat absorbed from the initially hotter sys- 
tem will be 



Q x = C(Tx - T a ). 



(27) 



Now we consider the situation in which we have lack of 
information about the exact values of these intermedi- 
ate temperatures T a and Tj,. We want to estimate the 
properties of this engine, taking into account the prior 
information we have about the parameters. Now it is 
clear that we have to assign the prior either to T a or 
Th, because the two parameters are related by Eq. (|25l) . 
Imagine two observers A and B who respectively take 
T a and Tb, as the uncertain parameter for the considered 
thermodynamic process. Thus from the perspective of A, 
the work is given by 



W = C[T 1 +T 2 -T a - 



while from £?'s point of view 



W = C[T 1 +T 2 -T b - 



TiT 2 



T{T 2 



(28) 



(29) 



Further we assume that A and B assign the same func- 
tional form for their priors and the range over which the 
parameters can take values is also the same. Thus prob- 
ability distribution P for the temperatures T a and Tf, is 
defined in the range [Ti,T2]. Now due to the constraint 
relating T a and T bl the probabilities assigned to any pair 
of values related by Eq. ([25]) . should be the same. This 
means 



P(T a )dT a = P(T b )dT b . 
Using (TJSD and we get: 

P(T a ) ' ' 



ln(l/60T a ' 



(30) 



(31) 



where = T 2 jT x . 

A. Estimated Work and Efficiency 

The work estimated by A is defined as: 



(W) 



To 



where W is given by Eq. 



(W) = CTi 



WP(T a )dT a , 

Thus 
2(l-i 



(1 + 6) 



In 9 



(32) 



(33) 



Similarly, the estimate for the heat absorbed is given by 

(1-9) 



(Qi) = CT X 



1 



In 9 



(34) 
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The efficiency, 77 = (W)/(Qi), is given by: 
01n0 + fl - 0) 



77 = 1 



ln0 



1 



(35) 



The observer B aiso arrives at the same estimates for 
work and efficiency. Now the expected values of interme- 
diate temperatures calculated by A and B after the work 
extraction process are 



(T a ) = (T b ) = Ti 



In (1/0)' 



(36) 



Moreover, all these estimates are also the same as derived 
from a quantum heat engine in the asymptotic limit in 
section IIVBI 



VI. CONCLUSIONS 

We have analysed the case of heat engines by assum- 
ing uncertainty in the exact values of the internal energy 
scales of the working medium. We have suggested the 
appropriate prior distributions for the uncertain param- 
eters based on prior information. In the case of quantum 
model, where the level spacings are the unknown param- 
eters, the expected values of work, heat and efficiency 
are equivalent to the expected values calculated from the 
classical model, where the intermediate temperature is 
the unknown parameter. Moreover in both the models, 
expected temperatures are equal at the end of work ex- 
traction. So the quantum model with two unknown pa- 
rameters shows similar behavior as the classical model, 
with a single uncertain parameter. In this point, it is in- 
teresting to analyse a special case of the quantum model 
with a single unknown parameter. Such kind of situation 
arises when the efficiency (77 = 1 — 02/01) of the engine 
is given. This simplifies the problem to that of a single 
uncertain parameter, either a\ or 02- Introducing two 
observers as discussed in Q, we get the functional form 
of the prior as 7r(a,) oc 1/a,. We can then calculate the 
expected work per cycle. In the asymptotic limit, the 
expression for expected work reduces to 



(37) 



Let us compare the above expression with its classical 
analog. The classical model of the heat engine discussed 
in section W\ has an efficiency 77 = 1 — T a /T\ = 1 — T2/T;,, 
when T a and T b are specified. Substituting this efficiency 
in Eqs. (|2"51) and (j2"9")l . we get the expression for work 
same as the quantum case given in Eq. (|37[) . 

The efficiency at the maximum expected work 
(Eq. ([37| ) is equal to 1 — y/0, the well known Curzon- 
Ahlborn efficiency [14] . In near equilibrium limit, this 
efficiency can be expanded as 



+ 0(1 -ef 



(38) 



So this efficiency appears at the maximum value of the 
expected work in the asymptotic limit for the quantum 
model with one uncertain parameter as well as at maxi- 
mum work for the classical model without any uncertain 
parameter. This efficiency falls in a certain universality 
class [lU [l(| , where the leading term in the near equilib- 
rium expansion is half of Carnot efficiency (r) c ). At this 
point, it is interesting to analyse the near equilibrium 
expansion of the efficiency expressed in Eqs. (1221) and 
(El): 



V 



(1-0) , (I- 
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[i -of 

135 



+ O[(l-0) 4 ]. (39) 



This expected efficiency is obtained in two cases, from 
the quantum model where two parameters (<ii and 02) 
are uncertain and also from the classical model with a 
single uncertain parameter (T a or T b ). The leading term 
in the near equilibrium expansion of this efficiency is 77c /3 
instead of 77 c /2 observed in Eq. (|3"51) . As another exam- 
ple, the efficiency at maximum power for an irreversible 
Brownian heat engine [I?} when optimized over the load 
and barrier height is given by 



V = 



2(1 



3-20(l + ln0) -0 2 ' 



(40) 



Expanding this efficiency for the near-equilibrium regime, 
we get 

,.,<l_^) + <i_^ + <i_J2! + 0|(1 _, n . (41) 

The similarity between Eqs. (|3"9")l and (|4"TT) up to the 



£ 0.4- 




FIG. 1: The efficiency in Eq. (|35|) and (l40[l are plotted versus 
9. The lower curve (thick line) is the efficiency at maximum 
power for the Brownian heat engine. The upper curve is the 
efficiency obtained from the classical model with intermediate 
temperature as the single uncertain parameter. The same 
expected efficiency is also obtained in the asymptotic limit of 
the quantum model when two intrinsic energy scales of the 
working medium are uncertain. 

second order, suggests there might be a new universality 
class which includes these efficiencies. In Fig. 1, the effi- 
ciencies given in (|3"5j) and P0|) are plotted. To conclude, 
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when all the unknown parameters are treated with prior 
probabilities, the expected behavior of quantum model 
with two internal energy scales as the uncertain param- 
eters is similar to the classical model with intermediate 
temperature as the single uncertain parameter. Interest- 
ingly, the expected efficiency in both cases, approaches 
?7 c /3 value, in the near-equilibrium limit. Further, the 
analysis based on subjective treatment of the incomplete 
information, suggests a new line of enquiry of a possible 
connection with optimal performance of the finite time 
irreversible models of heat engines. 
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